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Abstract 

Small- world networks (SWN), obtained by randomly adding to a regular struc- 
ture additional links (AL), are of current interest. In this article we explore (based 
on physical models) a new variant of SWN, in which the probability of realizing an 
AL depends on the chemical distance between the connected sites. We assume a 
power-law probability distribution and study random walkers on the network, fo- 
cussing especially on their probability of being at the origin. We connect the results 
to Levy Flights, which follow from a mean field variant of our model. 

PACS numbers: 05.40Fb, 05.60.-k, 71.55.Jv. 

1 Introduction 

Recently a lot of interest has centered on the so-called small-world networks (SWN) |l],U 
where an underlying regular lattice is supplemented with additional links (bonds), a fact 
which drastically reduces the minimal distances (the so-called chemical distances in the 
fractal literature) between pairs of points on the lattice |I],[3|-|8|. This question is of utmost 
importance for general network structures, say Internet-links ||, and for the spreading 
of diseases |3, 10-13 1, topics which depend on the minimal distances. On the other hand 



other questions are envisageable over such structures, for example random transport [1 I 



this requires solving diffusion-type problems, which are mathematically described by the 
Laplacian on the structure [I6j and the corresponding eigenvalues and eigenvectors |T] 



Examples of such problems are anomalous transport of charges and of excitations over 
networks |1§ . Most recent S WN-studies center on a one-dimensional chain supplemented 



with additional links (AL), which connect sites that are arbitrarily far from each other 
on the underlying lattice. While this being the simplest SWN envisageable, there are 
situations in which links between distant sites occur naturally; however, their lengths 
are then not necessarily uniformly distributed: Considering a polymer chain in solution, 
monomers which are far apart along the backbone can be quite close to each other in real 
space, so that for instance energy transfer over the structure may take crosscuts along 
sites near to each other in space |T8[ . Now the probability of having such close monomer 



pairs is related to the return to the origin of random walks, possibly under self-avoiding 
constraints. In this case the probability p(l) that two sites far apart along the backbone 
get close together in space, goes approximately as an inverse power-law of their mutual 



chemical distance I 19 



P(l) = £ (1) 

In this work, we will focus on networks constructed as in the SWN case, while however 
letting the additional bonds be distributed according to Eq. ([I]). We call these structures 
generalized small- world networks (GSWN). Clearly, the original SWN is recovered from 
the GSWN by letting a — > 0. On the other hand, in GSWN with a ^> 1 practically only 
sites which are already very close on the chain get to be connected; such GSWN have 
(apart from disorder) properties similar to the underlying regular lattice. Most interesting 
are GSWN with < a < 3, which show a wealth of features, because of the long-range 
character of the additional links. 

In this paper we will study random walks over GSWN, and especially the probability 
Po(t) of the walker being (still or again) at the origin; as discussed in previous works 
fT3], [2lf , this quantity reflects many of the properties of the density of eigenvalues of the 
underlying structure, and is easily obtainable by very effective, easy to program numerical 
procedures. 

Our paper is structured as follows: In the next section (Sec. ||) we discuss the con- 
struction of GSWN in more detail. In Sec. |3| we study the behavior of random walkers on 
GSWN. We find that for a well below 2 we have a behavior qualitatively similar to that 
of walkers over the SWN. However, for a larger than 2 we move towards another regime, 
quite reminiscent of random walks on regular lattices. The transition appears to happen 
around asi2, which prompts us to consider transient versus recurrent walks in Sec. |J In 



this Section we determine analytically Po(t) for a mean field variant of the GSWN model, 
which we then compare with the numerical findings Po(t) on GSWN. It turns out that 
the mean-field approach is related to Levy Flights and Levy Walks. Finally we close our 
paper by summarizing our conclusions in Sec.[5]. 

2 Construction of GSWN 

The construction of GSWN follows at first the SWN procedure closely: We start from a 
ring of N sites (i.e. a closed, regular Id lattice). Then we consider each site consecutively, 
and let it sprout with probability q an additional bond, which connects it to another site, 
see Fig. [TJ. We now let p(l) in Eq. (|I|) be the conditional probability that this bond gets 
attached to a site at the (minimal) chemical distance I from the sprouting site, measured 
along the ring. Here, given our periodic boundary conditions, the chemical distance I lies 
between 1 and int(iV/2), where int(x) denotes the largest integer X such that X < x (see 
Fig. |l|). Note that through p(l) our model differs from the standard SWN, where no I 
dependence is accounted for. In Eq. ([!]), the constant a normalizes the distribution so 
that 

int(N/2) 

2 J2 al ' a = l ( 2 ) 

i=i 

Now the exponent a in Eq. (P is a parameter of the model and will be varied in the 
following. For a finite system we can choose it freely, so that a > 0. In fact the choice 
a = 0, i.e. p(l) = 1/(N — 1) recovers one of the basic constructions of the SWN, by which 
connections to all sites but the source are equiprobable. For an infinite network, on the 
other hand, care has to be exercised; in order to keep Eq. ([!]) normalizable one has to 
have a > 1. 

We now turn to the basic procedure, in which the structure of the model enters through 
its connectivity matrix A. Now A is defined as follows: The off-diagonal elements of A, 
namely A^ with i ^ j, equal minus the number of links between the sites i and j of 
the network. The diagonal elements An count the total number of bonds connected to 
i. Hence the connectivity matrix is symmetric and one has Yli A%j = 0- Furthermore 
det(A) = and exactly one eigenvalue of A, say Ei, vanishes. One should note that 
from the spectrum of the A matrix one can determine e.g. the diffusion and vibrational 



properties of the structure, as well as its behavior in external fields [17,^1]. The spectrum 



of A for the SWN (a = 0) has been recently studied by Monasson |16j ; among his findings 
was the existence of a pseudo-gap in the density of states, a property which affects the 



long time diffusion properties PTJ. We note that A can be viewed as arising from two 



sources: One term, A^ , is deterministic and is due to the underlying regular lattice (here 
the ring). Another one, A^, is stochastic and arises from the randomly added links. Thus 



A = A^ + A^ 2) . Now formally A]^ = — 1, A^ = 2 and Ay = otherwise, where we 

identify site N + 1 with 1 (periodic boundary conditions). The entries in A\- are random, 
and for i ^ j equal 0, —1 or —2. In fact, letting I be the chemical distance between i and 
j {i ¥" 3)j one nas f° r the probability Pi(c) that A\, = —c 

Pi{c) = Q {qar a Y (1 - qar a ) 2 - c (3) 

The diagonal elements are, as before, determined from the requirement thatj^ A\- = 0. 
We close on a small note of caution by remarking that, due to our prescription, even for 
j G {1,... , i — 1}, the elements A\-' are not independent of each other. Thus if A\- = — 2 
for i ^ j, then for k ^ {i, j} one cannot have A ik = —2; a nondiagonal element having a 
value of —2 implies namely both for j and for k, that one of their additional bonds has 
started at i. By construction, however, i can only be the source of one additional bond. 
For decreasing q and increasing N, however, we expect such correlations between the A\, 
to be less and less important. 

3 Probability of Being at the Origin 

As a simple dynamical problem on the underlying structure, we focus now on the proba- 
bility for a random walker to be (still or again) at the origin of its walk at a later time. 
This quantity is fundamental for fractal lattices, where it leads directly to the spectral 
(harmonic) dimension [|15 |, a quantity of much importance ||22|| . As we have shown in a 



previous work, determining this quantity through a numerical cellular automaton proce- 



dure is quite straightforward and very revealing for SWN |14| . We look at the probability 
P(i,t\m) for the walker to be at site i at time t, given that it started at site m at time 
t = 0. One notes first that P(i,t\m) obeys the following master equation: 

»!) =- a ± Ali P(j, t \m), (4) 

j = l 

where a is a transition rate. In vector notation p( m )(£) = (P(l, t\m), . . . , P(N, t\m)) this 
relation reads: 

a -^i = -aA^\t). (5) 



and has the formal solution: 

P (m) (t) = exp (-a At) P (m) (0). (6) 

Now the initial condition is p( m )(0) = (0, . . . , 1, . . . ,0) with a single non-zero element at 
m. The probability that the walker is again at m at time t reads 

P(m,t\m) = \PM{t)] m = J>xp (-aAi)]^. [P^(0)],- = [exp (-<jAt)] mm (7) 

3 

This expression simplifies by averaging over all starting points, since then 

N _ , AT 



— Y^ p ( m > *M = ^ Tr ( ex P (-crAt)) = — J] e" 

m=l i=l 



.E;<rt 



holds, where Tr denotes the trace operation and ^ with 1 < % < N are the eigenvalues 
of the (symmetric) connectivity matrix A. Note that in Eq. @ because of the averaging 
over all initial points only the eigenvalues enter. Furthermore one can now readily average 
over different realizations, obtaining (since E\ = 0): 

N 

Po(t) = i + ^(E 6 "^^) = 4 + f P(E)e- Et dE (9) 

with p(-E) being the spectral density for £" > 0. 

We turn now to our calculations, by which we determine numerically Po(t) for different 
choices of a and q. For systems of size N = 1001 we construct the connectivity matrix, 
diagonalize it, and employ Eq. flD to evaluate Po(t). We use for each choice of a and q 
100 realizations to average over the structural disorder. In Fig. (fj) we display on double 
logarithmic scales Po(t) versus the dimensionless time at for q = 0.05 and for a ranging 
from to 3. First we note that for very long times, Po(t) reaches the constant value 1/N, 
which arises due to the eigenvalue E\ — 0. Increasing the size of the small- world network 
(i.e. N) pushes the long time plateau to lower values, but, as we have shown in an earlier 



work for a = [14|, leaves the Po(t) curves above the plateau practically unaffected. This 
is also what we find here for general a; this allows us to infer the qualitative features of 
Po(t) in the limit N — >• oo. 

Turning now to the discussion of the results, we note first that for a = they agree 
perfectly with our previous small- world network analysis ||14|| , which was based not on 



the diagonalization of A but on a cellular automaton method. The decay of Po(t) fol- 
lows at early times a power-law, which turns at later times into a stretched exponen- 
tial behavior. Asymptotically, the decay obeys to leading order the form exp(— Ct 1 ^ 3 ) 



H|.pT] (with C a constant), which follows from the spectral density of Ref. |TB| , p(E) ~ 
£ ,_1 ^exp(— O E~ x l 2 \ Note that in Fig. [| the curves flatten with increasing a, a sign 
that with growing a, a walker is less prone to get far away from its starting site. At early 
times Po(t) is little affected by variations in a, since at very short times it does not matter 
whether the AL bring the walkers very far away or not. The transition to pure power-law 
behavior appears to happen roughly around a = 2. 

Moving on to larger q to examine whether this transition depends on q, we plot in 
Fig. the results for q = 0.8. In this case there are more AL, and the results are more 
sensitive to the value of a: in Fig. |3] the cases a = 3.0 and a = 5.0 are easily distinguished. 
However, increasing a further does not change the curves significantly. The curves are 
also more spread out in the short-time domain than in the case q = 0.05. This is due to 
the fact that the quasi 1-dimensional behavior of the walk is mainly felt on distances of 
the order of 1/q, this being a measure of the mean distance between branching points fl4 |. 



We hasten to note that for very small times (not displayed on the Fig. ^) the curves for 
different a do coincide. Despite these differences, the qualitative behavior of the curves on 
Fig. H] and Fig. § is comparable. Furthermore, the cross-over behavior of a = 2 (given as a 
dotted line) appears even more clearly on Fig. |3|: the curves with a > 2 follow power-law 
decays closely, while the curves for a < 2 are partly concave, thus displaying a faster than 
power-law decay. 

4 Typology of Random Walks 

Let us briefly recall some terminology from the theory of random walks. A random walk 
is said to be recurrent, if it returns with probability 1 to the origin at some later time. 
Otherwise the walk is called transient. For a walk to be transient requires an infinite 
system, because in finite systems all walks (disregarding situations with traps, mortal 
walkers, etc) are recurrent. On homogenous lattices, a walk is transient if and only if 



dtP (t) (10) 

is finite. In line with this, we could expect walks on the GSWN with small a to be 
transient, given that for a stretched exponential behavior (j3 > 0): 



/oo 
dt exp (-Ct 13 ) < oo. (11) 



Moreover, in the opposite limit of large a, we observe that Po(t) follows a power-law decay 
with the exponent being nearly —1/2. It follows that for large a we have I = oo, an 



indication that the walk is recurrent. As we discuss in detail in the following, for walks on 
an infinite regular linear chain whose steps are long-ranged and obey Eq. (P for the step- 
lengths, the transition between recurrence and transience occurs at a = 2 P^| . It is now 
tempting to aim at explaining our findings of Sec. |3] along such lines. Such a connection 
is achieved by replacing the random small-world network-structure under investigation 
here by a regular one (a mean-field type approach), and letting the transition rates reflect 
the underlying statistics of the links |23| . In this way the probability of taking a step of 
length / > 1 is proportional to l~ a . However, as we show in the following, this regularizing 
approach is not particularly succesful, since it does not describe Po(t) well for small and 
moderate a. 

We start now from the so-called the Riemann walks p4 |, which are symmetric random 



walks on the linear chain, where each step of the walk can extend over the length / with 
probability 

7r(0 ~ l~ a , a > 1 (12) 

Such walks are recurrent for a > 2 and transient for a < 2 ||24|| . Riemann- walks are 



examples from the more general class of Levy Flights and Levy Walks [25] . Turning now 
to the problem of averaging both the GSWN structures and the random- walks over them, 
we simply replace in Eq. (||) A^ 2 ^ (remember that A = A^ 1 ) + A^ 2 ^) by its average (A^ 2 ^) 
over all GSWN. For the averaging we may use Eq. (E3) and obtain 

(A^) = 2al- a q = c(l), (13) 

where I is the chemical distance between i and j. By doing this we have now as connectivity 
matrix A = A^ 1 ) + (A^), whose eigenvalues Ek are readily found; they read (for iV odd): 

(JV-l)/2 

E k = 2 - 2 cos(2vr£;/iV) + 2 ^ (1 - cos(2nkj/N))c(j), (14) 

where k — . . . N — 1. Now Eq. fll4|) can be used to determine numerically Po(t) via Eq. 

In Fig. H we compare for iV = 1001, q = 0.05 and a = 3, a = 1.5 and a = the results 
of the two approaches. For a = 3 the two methods lead to a nice agreement; it seems that 
for a around or larger than 3 the fluctuations due to the disorder play only a minor role. 
On the other hand, as exemplified by a = 1.5 and a = 0, for a below 2 the mean field 
approach leads to Po(£)-f° rms which are quite different from those obtained in Sec. [| 



5 Conclusion 

In this work we have studied a new variant of the small-world network (SWN) model 
which takes into account the fact that the probability of adding links can depend on 
the chemical distance between the connected sites. Exemplarily, here we have taken the 
probability distribution to be a power law (with exponent a) of the chemical distance, see 
Eq. ([l|). We have focused on random walks and especially on the probability Po{t) of a 
random walker to be at its origin. Depending on the value of a we have found qualitatively 
different behaviors. Specifically, we found clues indicating that in the infinite system limit 
random walks on GSWN may change from being transient to being recurrent, as a crosses 
the marginal value of 2 from below. Moreover, we have shown that our model is related to 
Levy flights and to Riemann walks. We also found that a simple mean-field regularization 
of the GSWN-problem gives poor results for small a. Overall, it follows that GSWN with 
a < 2 are objects whose dynamical properties differ significantly from those of regular 
lattices. 
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Figure 1: Illustration of the small- world network. In this example two additional links 
are added, with the corresponding distances given on the figure. 
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Figure 2: The probability of being at the origin, Po(t), for q = 0.05 and a = 
3, 2, 1.5, 1.25, 1, 0.5 and 0.0 from upper to lower right. The curves for a = 2 and a = 3 
are hardly distinguishable in the figure. 
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Figure 3: Same as Fig. ||, but for a choice of q = 0.8. From upper to lower right are 
a = 5.0, 3.0, 2.2, 2.1, 2.0, 1.9, 1.8, 1.7, 1.5, 1.0, 0.5 and 0.0; the curve for a = 2 is indicated 
by dots. 
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Figure 4: Comparison of the mean field theory prediction (full curves) with the numerical 
data (dotted curves) for q = 0.05. The values of a are from above a = 3, a = 1.5 and 
a = 0. Good agreement is found only for large a, here a = 3. 
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